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Abstract 
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In  order  to  understand  the  physics  of  hydrogen  embrittlement  occurring  in 
Ti]_xAlx  alloys,  we  calculated  the  equilibrium  structure  and  the  elastic  con¬ 
stants  of  bulk  Ti,  Al,  the  ordered  alloy  T^Al,  as  well  as  their  hydrides,  using 
the  ab  initio  Density  Functional  formalism.  Calculations  for  ordered  systems 
were  performed  using  a  bulk  method  with  local  orbitals  and  ab  initio  pseu¬ 
dopotentials.  Results  for  the  electronic  structure  of  Ti  near  atomic  vacancies, 
substitutional  Al  and  interstitial  H  impurities  were  obtained  using  the  cluster 
approach  and  the  Linear  Muffin- Tin  Orbitals  in  the  Atomic  Spheres  Approxi¬ 
mation  (LMTO-ASA)-Green  function  formalism.  The  complementary  results 
of  these  calculations  provide  an  improved  understanding  of  the  problem,  why 
fracture  is  facilitated  upon  hydrogen  loading  in  Ti  crystals  containing  Al  im¬ 
purities. 
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I.  INTRODUCTION 


Titanium  is  a  very  important  metal  in  the  industry  due  to  the  diversity  of  its  usage 
in  marine,  aerospace  and  medical  implant  technology.  The  wide  range  of  its  applicability 
arises  chiefly  from  the  very  large  strength-to- weight  ratio  of  Ti.  It  is  also  well  known  that  the 
strength  and  other  mechanical  properties  of  this  system  improve  significantly  upon  alloying 
titanium  with  aluminum,  in  a  fashion  comparable  to  the  improved  elastic  behavior  of  steel 
caused  by  the  presence  of  carbon  in  the  iron  matrix.  Unfortunately,  the  ductility  of  the 
titanium  aluminides  TijAl  and  TiAl  is  very  low  below  600° C.  Moreover,  the  mechanical 
stability  of  Tii_xAlx  alloys  is  significantly  reduced  in  the  presence  of  hydrogen  [1]. 

Many  experimental  and  theoretical  studies  have  been  performed  during  the  last  twenty 
years  (for  a  review  see  Ref.  [2])  with  the  objective  to  reduce  the  brittle  behavior  of  these 
systems  at  room  temperature,  without  deteriorating  their  advantageous  elastic  behavior  at 
elevated  temperatures.  These  studies  have  established  that  the  origin  of  the  low  ductility 
at  room  temperature  lies  in  the  hindrance  of  motion  of  specific  dislocations  due  to  their 
pinning  [3],  and  that  this  effect  is  related  to  the  flow  stress  anomaly  found  in  TiAl  single 
crystals  [4]. 

Two  theoretical  calculations  have  addressed  the  pinning  of  dislocations  in  the  stoichio¬ 
metric  TiAl  alloys,  both  of  which  relate  this  effect  to  the  covalent  bonding  in  these  systems. 
Greenberg  and  co-workers  have  shown  that  the  anisotropy  of  the  charge  density  near  the 
Ti  atoms  in  TiAl  causes  a  many-valleyed  Peierls  relief  which  could  affect  the  mobility  of 
dislocations  [5].  More  recently,  Morinaga  and  coworkers  found  a  strong  Ti3d-A13p  bonding 
between  the  (001)  planes  of  the  TiAl  alloy  [6],  which  may  lie  at  the  origin  of  the  large  value 
of  the  sheax  modulus  C44,  determined  by  Fu  and  Yoo  [7].  Further  alloying  with  3 d  transi¬ 
tion  metals  might  reduce  this  strong  Ti3d-A13p  bonding,  thereby  enhance  the  mobility  of 
dislocations  and  hence  the  ductility  of  TiAl. 

The  objective  of  our  study  is  to  understand  the  mechanical  properties  of  the  more  dilute 
Tii_xAlx  systems  and  the  corresponding  alloys  containing  hydrogen  [8].  Even  at  low  A1 
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concentrations,  the  hydrogen-free  Tii_xAlx  alloys  show  an  increased  strength  as  compared 
to  pure  Ti.  Unfortunately,  this  increased  strength  is  accompanied  by  an  enhanced  propensity 
towards  stress-corrosion  cracking  (SCC)  in  aqueous  solutions  [1,9]  as  well  as  organic  liquids 
[10]. 

Our  discussion  will  focus  on  the  transgranular  (rather  than  the  intergranular)  cracking 
mechanism  which  dominates  especially  at  larger  A1  concentrations  used  in  high-strength 
alloys.  Even  though  the  onset  of  transgranular  fracture  can  be  influenced  by  many  factors, 
the  reduced  stress  resistance  has  been  traced  back  in  many  cases  to  hydrogen  uptake  by 
the  metal.  Presence  of  hydrogen  has  been  found  to  cause  significant  strength  reductions 
in  Ti  alloys  containing  as  little  as  5%  (weight)  of  A1  [9].  It  is  interesting  to  note  that 
the  pure  Ti  metal,  which  forms  a  stable  hydride,  is  not  susceptible  to  hydrogen-induced 
transgranular  stress-corrosion  cracking,  and  neither  are  Ti-Al  alloys  containing  less  than  2% 
(weight)  Al.  It  has  also  been  observed  that  with  increasing  A1  concentration  in  Ti,  the  nature 
of  dislocation  motion  changes  from  “tangled”  or  cellular  to  a  co-planar  slip  at  essentially 
the  same  composition  which  marks  the  onset  of  the  susceptibility  towards  transgranular 
stress-corrosion  cracking  [11]. 

These  observations  raise  a  number  of  issues  which  have  to  be  addressed.  The  first  basic 
question  is,  whether  the  adverse  effects  of  hydrogen  on  the  mechanical  stability  of  Tii_xAlx 
alloys  with  x  >  2%  has  a  predominantly  electronic  origin,  reflected  in  the  change  of  local 
bonding  (in  a  way  similar  to  the  effect  of  H  on  the  stability  of  Pd  [12]),  or  rather  a  structural 
origin,  related  to  a  propensity  towards  structural  transitions. 

An  indication  in  favor  of  the  electronic  origin  are  the  significantly  different  electronic 
configurations  and  electronegativities  of  Al  and  Ti,  implying  a  large  charge  transfer  between 
these  species  in  the  alloy.  The  interrelationship  between  the  charges  carried  by  the  atoms  and 
mechanical  failure  has  been  demonstrated  by  a  successful  application  of  cathodic  potentials 
to  inhibit  transgranular  cracking  [13]. 

On  the  other  hand,  Tii_xAlx  is  not  a  random  alloy  even  at  the  low  Al  concentrations 
associated  with  the  onset  of  hydrogen  embrittlement;  at  higher  Al  concentrations  and  tem- 
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peratures  the  propensity  towards  the  formation  of  the  intermetallic  compound  T^Al  has 
been  observed.  This,  together  with  the  fact  that  the  crystalline  structures  of  the  pure  com¬ 
ponents  are  different,  would  rather  speak  in  favor  of  a  structural  origin  of  stress-corrosion 
cracking.  Hence  one  of  the  first  questions  to  be  asked  is  whether  it  is  predominantly  elec¬ 
tronic  or  structural  effects  which  lie  at  the  root  of  hydrogen-induced  stress- corrosion  cracking 
in  these  systems. 

To  our  knowledge,  no  theoretical  information  is  available  on  the  electronic  and  structural 
properties  of  Tii_xAlxHy  alloys,  whereas  several  studies  have  been  performed  on  Ti  hydrides 
[14].  We  have  made  use  of  three  different  yet  complementary  ab  initio  techniques  in  order 
to  elucidate  different  aspects  of  hydrogen  embrittlement  in  Tii_xAlx  alloys  at  various  com¬ 
positions.  The  objective  of  our  studies  is  not  to  compare  the  range  of  applicability  of  these 
techniques,  but  rather  to  obtain  a  detailed  and  reliable  description  of  the  bonding  changes 
in  the  Tii_xAlx  matrix  in  the  vicinity  of  hydrogen. 

A  summary  of  the  theoretical  methods,  which  have  been  described  in  more  detail  else¬ 
where,  will  be  presented  in  the  following  Section  II.  Results  regarding  the  pure  metals  A1 
and  Ti  will  be  given  in  Section  III.  In  Section  IV,  we  present  results  for  the  hydrogen-free 
Tii_xAlx  alloys,  and  in  Section  V  for  these  alloys  with  dissolved  hydrogen.  We  discuss  our 
results  in  Section  VI  and  present  our  conclusions  in  Section  VII. 

II.  THEORY 

We  determine  the  equilibrium  geometry,  the  total  energy,  and  the  electronic  structure  of 
hydrogen-free  and  hydrogen-loaded  Tii_xAlx  alloys  using  different  implementations  of  the 
Density  Functional  formalism  within  the  Local  Density  Approximation  (LDA)  [15].  In  the 
electronic  ground  state,  the  total  energy  functional  is  minimized  by  solving  self-consistently 
the  set  of  Kohn-Sham  equations  [15] 

h2 

[—  „ -  V2  +Von  +  Vh  +  VrcjV’n  =  tn'&n  (l) 

1  m 

and 
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(2) 


pif)  =  Z)IV’njt(^)l2  • 

nk 

Here,  Vlon  is  the  potential  due  to  the  ion  cores,  for  which  we  use  either  the  all-electron  po¬ 
tential  or  an  ab  initio  pseudopotential.  V#  is  the  Hartree  potential  and  Vxc  is  the  exchange- 
correlation  potential  associated  with  the  electron  density  p(r),  which  we  use  in  the  form 

parametrized  by  Hedin  and  Lundqvist  [16].  and  en£  are  the  Kohn-Sham  wavefunctions 

— * 

and  energies,  respectively,  where  the  crystal  momentum  k  is  used  only  in  the  case  of  ordered 
alloys.  The  wave  functions  are  expanded  in  a  local  basis,  containing  atom-centered  and 
floating  orbitals,  and  the  angular  momentum  is  limited  to  lmax  =  2  in  all  our  calculations. 
The  total  energy  of  the  system  in  the  electronic  ground  state  is  calculated  using 


Etot  =  Yj  enk  ~  \  J  VH{pif))pif)  ~  J  V*c[p{r))p(r)di 

nk 

+  J  dr  exc[p(r)]p(r)  +  E 

ion— ion  5 


(3) 


where  e;rC[/j(f)]p(f)  is  the  exchange-correlation  energy  density  and  E ion-ion  describes  the 
interaction  between  ion  cores. 


A.  Calculations  for  ordered  bulk  alloys 

To  calculate  the  electronic  properties  of  ordered  bulk  alloys,  we  use  the  pseudopotential 
local  orbital  method  which  has  been  discussed  in  detail  previously  and  which  has  been 
applied  successfully  to  transition  metals  and  their  hydrides  [12,17].  The  interaction  of  valence 
electrons  with  the  ion  cores  is  described  by  first-principles  pseudopotentials  Vion  of  the 
Hamann-Schliiter-Chiang  type  [18].  The  Gaussian  local  orbital  basis  functions 

faim( f)  =  Aalme-ar2Klm(6,  <j> ),  (4) 

are  centered  on  atomic  sites,  and  are  also  used  as  floating  orbitals  at  interstitial  sites.  The 
Kubic  harmonics  Kim  represent  the  angular  part  of  the  basis  functions,  and  Aa[m  are  nor¬ 
malization  constants.  The  Brillouin  zone,  which  is  defined  in  ordered  crystalline  solids,  is 
sampled  by  a  homogeneous  mesh  of  special  ^-points  [19]. 
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B.  Calculations  for  ordered  solids  with  impurities 


To  obtain  the  charge  distribution  around  an  impurity  dissolved  in  an  otherwise  unper¬ 
turbed  medium  we  use  the  Linear  Muffin  Tin  Orbitals  (LMTO)-Green  function  method 
within  the  Atomic  Sphere  Approximation  (ASA).  This  is  an  all-electron  technique  which 
is  used  within  the  frozen-core  approximation.  In  the  multiple  scattering  theory,  the  local 
density  of  states  at  each  location  r  in  the  alloy  is  given  by  the  imaginary  part  of  the  Green 
function  G(f,f,E).  It  has  been  well  established  how,  within  the  LMTO  scheme,  this  func¬ 
tion  can  be  obtained  from  the  structural  Green  function  T(£')  [20],  which  is  related  to  the 
Green  function  of  the  pure  host  r(E)  by  the  Dyson  equation 

TTu(E)  =  TTkE)  -  E  r TkE)  A  Pp:'{E)Tt%{E) .  (5) 

n",L" 

Here,  to  labels  the  site,  L  =  (l,m)  the  orbital  symmetry,  and  A Pp"(E)  describes  the  local 
perturbation  at  the  site  to"  due  to  the  presence  of  the  defect.  Sites  at  which  A Pfi'(E)  is 
nonzero  define  the  perturbed  cluster.  Due  to  the  good  screening  in  metals,  the  perturbation 
is  essentially  confined  within  the  nearest  neighbor  shell  of  the  impurity  site,  which  justifies 
the  selection  of  these  sites  as  the  perturbed  cluster  in  our  calculations.  We  perform  the 
calculations  using  the  Simple  Diagonal  Cluster  (SDC)  model  in  which  the  intersite  Green 

_ nnn 

function  matrix  elements  of  the  pure  host  TLL„(E)  (to  ^  to")  are  neglected  in  Eq.  (5).  The 
SDC  approximation,  while  reducing  the  computing  time  by  one  order  of  magnitude,  has 
been  shown  to  yield  results  which  are  very  close  to  those  obtained  by  an  exact  solution  of 
Eq.  (5)  [20]. 

In  our  calculations,  we  assume  that  all  atoms  have  a  spherical  muffin-tin  potential  within 
their  Wigner-Seitz  sphere,  and  that  the  radius  of  these  overlapping  spheres  is  the  same  at 
all  the  sites  in  the  alloy.  Since  the  screening,  which  occurs  outside  the  perturbed  cluster 
in  response  to  the  impurity,  does  not  change  significantly  the  charge  distributions  in  the 
perturbed  cluster,  we  ignored  this  screening  effect  in  our  calculation. 
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C.  Cluster  calculations 


In  order  to  closely  investigate  the  local  impact  of  hydrogen  on  the  metallic  bonding 
in  Ti-Al  systems  with  no  symmetry  restrictions,  we  used  a  cluster  implementation  of  the 
Density  Functional  formalism  in  the  DMol  code  [21].  The  molecular  orbitals  are  represented 
by  linear  combinations  of  numerically  generated  basis  functions.  In  the  double-numerical 
basis  set  [21],  the  occupied  valence  orbitals  are  represented  by  two  wave  functions  each.  The 
basis  is  expanded  at  each  site  using  polarization  functions  with  angular  momentum  quantum 
number  which  is  by  one  higher  than  the  highest  angular  momentum  in  the  free  atom.  We 
use  the  frozen  core  approximation  to  treat  the  inner  core  electrons  and  the  von  Barth  and 
Hedin  potential  to  describe  the  exchange- correlation  energy  [22]. 

III.  PURE  CONSTITUENTS 

A.  Aluminum 

Even  though  the  electronic  structure  of  A1  has  been  studied  previously  [23] ,  we  include 
our  corresponding  results  for  the  sake  of  completeness  and  comparison  with  Al-based  alloys. 

The  energy  of  the  A1  pseudoatom  amounts  to  E^^m(Al)  =  —65.562  eV,  which  includes 
the  calculated  spin-polarization  energy  of  —0.178  eV.  The  3s  and  3 p  valence  orbitals  are 
represented  by  three  radial  Gaussian  functions  each,  with  the  optimized  decays  (in  a.u.)  a  = 
0.10,0.28,0.80,  yielding  12  independent  basis  functions  at  each  site.  In  our  bulk  calculation, 
the  Fourier  expansion  of  the  charge  density  is  truncated  at  an  energy  of  49  Ry,  and  the 
irreducible  part  of  the  Brillouin  zone  is  sampled  by  60fe-points. 

The  calculated  properties  of  bulk  A1  in  the  face-centered  cubic  (fee)  phase,  which  axe 
summarized  in  Table  I,  are  in  excellent  agreement  with  experimental  data  [24].  We  find  bulk 
aluminum  to  be  slightly  overbound  [25],  due  to  the  well  known  underbinding  of  atoms  in 
LDA.  The  total  density  of  states  (TDOS)  [26]  of  Al,  shown  in  Fig.  1(a),  exhibits  the  nearly 
parabolic  behaviour  near  the  lower  valence  band  edge,  as  expected  for  a  near-free  electron 
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metal.  We  found  these  results  to  lie  very  close  to  our  parallel  calculations  using  a  plane- wave 
basis  and  the  all-electron  LMTO  method,  indicating  the  convergence  and  interchangeability 
of  these  techniques. 


B.  Titanium 

The  analogous  calculations  for  Ti  yield  the  value  E°f£m[Ti')  =  —89.963  eV  for  the  total 
energy  of  the  isolated  pseudoatom.  The  Ti  atom  contains  both  localized  3 d  and  delocalized 
4s  electrons,  which  can  only  be  described  by  relatively  deep  pseudopotentials.  In  our  calcu¬ 
lation,  the  3 d,  4s,  and  4p  valence  orbitals  are  represented  by  four  radial  functions  each,  with 
the  optimized  Gaussian  decays  (in  a.u.)  a  =  0.18,0.56,1.79,  and  5.70.  Hence,  our  basis 
contains  40  independent  functions  at  each  Ti  site. 

The  equilibrium  structure  of  Ti  at  low  temperatures  is  hexagonal  close-packed  (hep), 
which  differs  from  the  fee  structure  only  in  the  stacking  sequence  of  the  close-packed  layers 
(ABAB  in  hep  vs.  ABCABC  in  fee).  Since  all  coordination  numbers  and  hence  the  packing 
fractions  in  the  two  structures  are  the  same,  we  expect  the  total  energy  and  the  electronic 
structure  of  Ti  in  these  two  structures  to  be  very  similar.  This  is  indeed  confirmed  by  the 
results  of  our  calculations,  given  in  Table  I  and  Figs.  1(b)  and  (c).  The  gain  of  symmetry 
(and  consequently  reduction  of  the  computing  requirements)  in  the  cubic  structure  more 
than  compensate  for  the  small  error  introduced  when  considering  the  fee  in  lieu  of  the  hep 
phase. 

These  bulk  results  are  based  on  a  cutoff  of  the  Fourier  expansion  of  the  charge  density 

and  potential  at  an  energy  of  49  Ry,  and  a  sampling  of  the  irreducible  part  of  the  Brillouin 

— *  *»♦  _ 

zone  by  27  k  points  in  the  hep  phase  and  60  k  points  in  the  fee  phase.  Our  results  for  the 

electronic  structure  of  Ti  in  the  hep  structure  lie  close  to  the  LMTO  results,  obtained  using 

196  k  points  in  the  irreducible  wedge  of  the  Brillouin  zone,  in  particular  reproducing  the 

pronounced  depression  in  the  density  of  states  near  Ep  discussed  below. 

A  close  inspection  of  the  structural  data,  presented  in  Table  I,  confirms  the  high  precision 
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with  which  LDA  predicts  the  equilibrium  geometry  of  bulk  Ti.  As  suggested  above,  the 
binding  energy  of  Ti  in  the  hep  and  the  fee  phase  are  nearly  the  same  [25],  and  found 
to  be  indistinguishable  in  our  LDA  calculations.  Our  total  energies,  on  the  other  hand, 
underestimate  the  cohesion  of  Ti  in  comparison  with  the  experimental  data.  In  view  of  the 
computational  difficulties  traditionally  associated  with  the  nonlocal  3 d  pseudopotentials,  we 
checked  the  convergency  of  these  results  by  augmenting  the  basis  set  by  interstitial  floating 
orbitals,  and  by  extending  the  Fourier  cutoff  energy  from  49  Ry  to  100  Ry,  but  found  the 
total  energy  to  be  affected  by  less  than  0.1  eV.  More  important  to  us  is  the  elastic  response 
of  Ti  where,  in  the  case  of  the  bulk  modulus,  the  LDA  results  are  in  close  agreement  with 
the  experimental  data. 

The  small  difference  between  Ti  in  the  hep  and  the  fee  structure  is  reflected  in  the 
electronic  densities  of  states,  which  are  presented  in  Figs.  1(b)  and  (c),  respectively.  The 
width  of  the  occupied  band  Wocc  ~  3  eV,  the  total  3d-bandwidth  Wtot  ~  6  eV,  and  the 
three-peak  structure  of  the  3d  band  below  the  Fermi  level  are  common  to  both  hep  and  fee 
titanium.  The  slight  stabilization  of  the  hep  with  respect  to  the  fee  phase  of  Ti  is  reflected 
in  a  state  depletion  in  the  hep  structure  near  the  Fermi  level,  which  is  accompanied  by  a 
lowering  of  the  band  structure  energy. 

IV.  HYDROGEN-FREE  Tlx.xALx  ALLOYS 

Before  discussing  the  adverse  effect  played  by  hydrogen  on  the  elastic  behavior  of  high- 
strength  Tii_xAlx  alloys,  we  discuss  in  detail  the  electronic  structure  of  the  hydrogen-free 
ordered  alloy  TisAl,  and  compare  the  results  to  those  for  a  substitutional  A1  impurity  in 
Ti(hcp). 


A.  Ordered  alloy  TijAl 

We  discuss  the  ordered  compound  T^Al,  since  the  presence  of  this  phase  in  high  strength 
titanium  alloys  has  been  linked  to  their  increased  susceptibility  towards  transgranular  crack- 
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mg  [27].  The  close-packed  structure  of  Ti3Al  which  we  consider  has  an  underlying  simple 
cubic  (sc)  Bravais  lattice,  and  can  be  obtained  from  the  metastable  fee  phase  of  Ti  by  replac¬ 
ing  the  corner  sites  in  the  conventional  unit  cell  by  A1  atoms.  This  structure,  which  is  very 
closely  related  to  that  of  the  experimentally  observed  hep  stacked  system  [28],  has  signifi¬ 
cantly  higher  symmetry  and  hence  can  be  computationally  treated  with  a  higher  precision. 
Our  basis  contains  Ti-  and  Al-centered  orbitals,  as  discussed  in  the  preceding  Section.  We 
found  our  results  to  be  well  converged  when  sampling  the  irreducible  part  of  the  Brillouin 
zone  by  20  fc-points,  and  when  using  an  energy  cutoff  of  49  Ry  in  the  Fourier  expansion  of 
the  charge  density. 

The  negative  value  of  the  formation  energy  of  Ti3Al  [29],  reported  in  Table  II,  indicates 
stability  with  respect  to  decomposition  into  pure  Ti  and  A1  metals.  The  calculated  and 
observed  equilibrium  structural  properties,  as  well  as  the  bulk  modulus,  of  Ti3Al  are  sum¬ 
marized  in  Table  I  [25].  The  experimental  value  given  for  the  lattice  constant  a  of  Ti3Al  in 
the  metastable  cubic  phase  has  been  obtained  by  extrapolating  the  experimental  data  for 
the  hexagonal  phase  [30],  assuming  that  the  unit  cell  volume  is  the  same  in  the  cubic  and 
the  hexagonal  structure.  This  value  is  in  good  agreement  with  our  calculated  equilibrium 
lattice  constant. 

The  electronic  density  of  states  of  the  ordered  Ti3Al  alloy  [26],  calculated  using  ab  initio 
pseudopotentials  and  a  local  orbital  basis,  is  presented  in  Fig.  2(a).  The  site-projected  local 
electronic  densities  of  states  (LDOS)  on  the  A1  and  Ti  sites  are  given  in  Figs.  2(b)  and  (c), 
respectively.  When  comparing  this  density  of  states  to  that  of  the  pure  bulk  constituents 
A1  and  Ti  in  Fig.  1,  we  notice  the  formation  of  a  new  hybrid  band  located  ~  5  eV  below 
the  Fermi  level,  with  a  bandwidth  of  about  1  eV.  The  angular  momentum  decomposition 
shows  that  this  band  has  predominantly  an  Al3s-Ti3d  character;  a  comparison  between  the 
local  densities  of  states  in  Figs.  2(b)  and  (c)  indicates  a  predominant  localization  of  these 
states  on  the  A1  sites.  The  presence  of  new  occupied  Al-derived  states  is  consistent  with 
the  Mulliken  population  analysis,  reported  in  Table  III,  which  indicates  a  net  transfer  of  as 
much  as  0.9  electrons  from  Ti  to  the  A1  sites. 
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B.  A1  impurities  in  Ti 


In  order  to  understand  the  effect  of  aluminum  on  the  electronic  structure  of  Tii_xAlx 
alloys  even  at  very  low  A1  concentrations,  we  performed  calculations  for  a  substitutional  A1 
impurity  in  a  pure  Ti(hcp)  matrix.  We  placed  all  the  atoms  on  a  perfect  Ti(hcp)  lattice, 
given  by  the  experimentally  observed  lattice  parameters  (a  =  5.577  a.u.  and  c/a  =  1.588). 
In  the  LMTO-ASA  calculation,  we  assumed  the  same  mean  radius  for  all  the  overlapping 
atomic  Wigner-Seitz  spheres,  and  considered  the  perturbed  cluster,  consisting  of  the  A1 
impurity  and  its  twelve  nearest  neighbors  [see  Fig.  3(a)],  to  be  embedded  in  an  otherwise 
perfect  crystal.  These  results  were  compared  to  those  obtained  applying  the  DMol  approach 
to  the  isolated  Ti^Al  core  cluster,  frozen  in  the  bulk  geometry. 

A  three-dimensional  representation  of  the  total  charge  density  within  the  isolated  Ti^Al 
core  cluster  is  shown  in  Fig.  3(b).  A  two-dimensional  contour  plot  of  the  charge  density 
in  the  (0001)  plane  containing  the  A1  impurity,  indicated  by  a  dashed  line  in  Fig.  3(a),  is 
presented  in  Fig.  3(c)  for  the  isolated  cluster  and  in  Fig.  3(d)  for  the  Ti  matrix  with  a  single 
A1  impurity.  A  comparison  of  the  results  in  Figs.  3(c)  and  (d)  indicates  that  the  isolated 
cluster  gives  a  good  idea  of  the  electronic  structure  near  a  substitutional  A1  impurity  in  bulk 
Ti. 

The  local  densities  of  states  on  the  A1  impurity  and  on  its  nearest  Ti  neighbors,  obtained 
using  the  LMTO-ASA  approach,  are  shown  in  Figs.  4(a)  and  (b),  respectively.  Our  analysis 
of  the  net  charges  on  the  A1  impurity  and  its  Ti  neighbors,  presented  in  Table  III,  yields  a 
consistent  picture  with  results  obtained  in  the  ordered  Ti3Al  alloy.  As  seen  in  Fig.  4(a),  the 
substitutional  A1  impurity  induces  a  sharp  resonant  state  near  the  bottom  of  the  valence 
band,  about  7  eV  below  the  Fermi  level.  This  state  is  mainly  localized  on  the  A1  site  and 
nearly  invisible  at  the  neighboring  Ti  sites  and,  in  the  ordered  Ti3Al  alloy,  spreads  out  to  the 
A13$-Ti3d  derived  band  discussed  above.  As  in  the  ordered  alloy,  the  A1  site  is  negatively 
charged.  The  LMTO-ASA  calculation  yields  for  the  total  charge  on  the  A1  site  the  value 
of  0.18  electrons,  which  compares  well  with  the  Mulliken  charge  of  0.22  electrons  on  the  A1 
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site  in  the  isolated  Ti^Al  cluster,  yet  is  significantly  smaller  than  the  value  in  the  ordered 
Ti3Al  alloy. 

V.  HYDROGEN-LOADED  Tlx.xALx  ALLOYS 

In  an  attempt  to  enhance  our  understanding  of  the  microscopic  origin  of  hydrogen  em¬ 
brittlement  in  Tii_xAlx  alloys,  we  calculated  the  structural  and  electronic  properties  of  bulk 
Al,  Ti  and  Ti3Al  at  different  hydrogen  concentrations,  with  the  hydrogen  atoms  located  in 
the  interstitial  sites  of  the  lattice. 

The  basis  functions  on  the  metal  atom  sites  have  been  described  above.  Hydrogen  atoms 
in  a  crystalline  environment  are  represented  by  sixteen  independent  basis  functions.  The 
radial  dependence  of  the  basis  functions  with  s  and  p  symmetry  is  given  by  four  Gaussians 

each,  with  the  decay  constants  (in  a.u.)  a  —  0.20,0.42,0.68  and  1.36  [12].  We  use  the  same 

«*♦ 

k  point  sampling  of  the  irreducible  Brillouin  zone  and  the  same  energy  cutoff  of  the  Fourier 
expansion  of  the  charge  density  in  hydrogen-loaded  as  in  the  hydrogen-free  systems. 

A.  A1H 

It  is  well  known  experimentally  that  hydrogen  is  not  readily  absorbed  in  metallic  alu¬ 
minum  [8].  For  the  sake  of  completeness  of  our  study,  nevertheless,  and  in  order  to  demon¬ 
strate  the  adequacy  of  our  computational  scheme,  we  present  in  the  following  results  for  a 
hypothetical  compound  A1H,  where  the  H  atoms  occupy  the  octahedral  interstitial  sites  in 
the  fee  lattice  of  aluminum. 

The  positive  calculated  value  for  the  formation  energy  of  A1H  [31],  given  in  Table  II, 
agrees  indeed  with  the  experimental  finding  that  A1H  should  be  unstable.  Our  predictions 
for  the  other  structural  properties  of  A1H  are  reported  in  Table  I.  These  results  suggest  that 
the  bulk  modulus  of  A1H  should  be  nearly  the  same  as  that  of  pure  Al,  whereas  the  shear 
modulus  C44  of  Al  would  be  expected  to  decrease  by  more  than  50%  in  the  hydride.  As 
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shown  in  Table  III,  this  change  of  the  elastic  response  is  accompanied  by  a  large  transfer  of 
0.7  electrons  from  the  A1  to  the  H  sites. 


B.  TiHx 

In  contrast  to  aluminum,  bulk  titanium  does  form  stable  hydrides  [8].  As  mentioned 
above,  we  constrain  our  structural  optimization  of  TiHx  to  the  fee  phase,  and  compare  our 
results  to  those  for  Ti(fcc).  Guided  by  available  information  [8],  we  place  the  hydrogen 
atoms  at  the  octahedral  and  tetrahedral  interstitial  sites  of  the  fee  lattice. 

We  first  consider  the  ordered  compound  TiH,  where  hydrogen  atoms  occupy  all  the 
available  octahedral  sites,  as  shown  in  Fig.  5(a).  The  observed  readiness  of  Ti  to  absorb 
hydrogen  is  reflected  in  a  negative  value  of  the  formation  energy  of  TiH  [31],  given  in 
Table  II.  The  structural  properties  of  TiH  are  summarized  in  Table  I.  Our  results  indicate 
an  expansion  of  6%  of  the  lattice  constant  of  Ti  upon  hydrogen  uptake,  in  general  agreement 
with  lattice  dilatations  observed  in  compact  metallic  structures  [8].  Results  of  our  Mulliken 
population  analysis,  reported  in  Table  III,  suggest  a  small  positive  charge  on  the  Ti  atoms. 
The  total  charge  density,  calculated  for  the  isolated  Ti^He  core  cluster  [see  Fig.  5(a)]  using 
the  DMol  formalism,  is  presented  in  Fig.  5(b). 

Next  we  consider  the  compound  TiH2,  where  all  the  hydrogen  atoms  occupy  the  tetrahe¬ 
dral  interstitial  sites  of  the  fee  lattice,  yielding  a  fluorite  structure  of  the  system.  As  shown 
in  Table  II,  the  calculated  formation  energy  [31]  of  TiH2  is  negative  and  lower  than  that 
of  TiH.  This  is  in  agreement  with  experimental  results  which  suggest  that  transition  metal 
elements  of  the  groups  3B,  4B,  and  partly  those  of  group  5B  tend  to  form  dihydrides  [32], 
The  structural  properties  of  TiH2,  summarized  in  Table  I,  suggest  a  further  expansion  of  the 
lattice  with  respect  to  TiH,  yet  essentially  the  same  compressibility  as  the  monohydride. 

We  find  that  our  calculated  total  density  of  states  N(E)  of  TiH2,  shown  in  Fig.  6(b),  lies 
very  close  to  Augmented  Plane  Wave  (APW)  results  for  the  same  system  [33].  In  particular, 
our  value  of  the  density  of  states  at  the  Fermi  level  N(Ep)  =  1.76  states/eV  [26]  compares 
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very  well  with  the  APW  value  N(Ef )  =  1.73  states/eV  [33].  A  comparison  between  the 
local  density  of  states  on  the  Ti  sites  in  TiH2,  shown  in  Fig.  6(c),  with  the  density  of  states  of 
pure  Ti,  reproduced  in  Fig.  6(a),  suggests  that  the  strong  Ti-H  interaction  modifies  strongly 
the  Ti3d  band  and  induces  a  new  Ti3d-Hls  derived  band  about  5  eV  below  Ep.  This  new 
H-induced  band  is  split  into  two  subbands  due  to  the  close  proximity  of  the  two  hydrogen 
atoms  in  the  unit  cell.  The  overscreening  of  the  proton  charge  on  the  H  sites,  reported  in 
Table  III,  is  much  smaller  than  in  TiH,  and  agrees  well  with  previous  results  [34]. 

C.  Ti!_xAlxHy 

In  order  to  better  understand  the  adverse  effect  of  hydrogen  on  the  binding  in  Tii_xAlx 
alloys,  we  performed  Density  Functional  calculations  for  the  ordered  compound  Ti3AlH4 
in  the  simple  cubic  structure.  We  compared  these  results  to  those  for  the  hydride  of  a 
very  dilute  TiAl  alloy,  represented  by  an  isolated  Tii2AlH6  core  cluster  shown  in  Fig.  7(a). 
Even  though  the  simple  cubic  structure  of  the  ordered  alloy  may  not  correspond  to  the 
equilibrium  phase  of  this  system,  we  selected  it  for  its  high  symmetry  and  for  the  sake 
of  simple  comparison  with  our  results  for  the  hydrogen-free  Ti3Al  compound.  As  in  TiH 
and  A1H,  we  assumed  that  the  hydrogen  atoms  occupy  all  the  octahedral  interstices  in  this 
structure,  four  in  number  in  the  conventional  unit  cell  of  the  fee  lattice  which  coincides  with 
the  unit  cell  of  Ti3AlH4.  We  expect  our  results  for  this  system  to  elucidate  the  conjugate 
roles  played  by  both  H  and  A1  in  the  embrittlement  of  Ti-based  systems. 

The  calculated  positive  value  of  the  formation  energy  of  Ti3AlH4,  given  in  Table  II, 
indicates  a  potential  instability  of  this  system  towards  decomposition  into  its  components 
in  their  standard  state.  This  may  be  the  reason  why  no  structure  has  been  observed  or 
reported  for  this  system  [35].  The  calculated  structural  properties,  summarized  in  Table  I, 
show  among  others  that  LDA  slightly  underestimates  the  experimental  equilibrium  lattice 
constant  [35].  The  compactness  of  this  structure  is  reflected  in  the  large  value  of  the  bulk 
modulus.  We  also  note  the  substantial  decrease  of  the  shear  modulus  C 44  of  this  system 
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with  respect  to  the  hydrogen-free  structure. 

The  calculated  total  density  of  states  [26]  of  bulk  Ti3AlH4  is  compared  to  that  of  the 
hydrogen-free  Ti3Al  system  in  Fig.  8.  The  corresponding  results  of  the  Mulliken  population 
analysis  for  these  bulk  compounds  are  reported  in  Table  III.  The  most  prominent  difference 
in  Ti3AlH4  with  respect  to  the  hydrogen-free  compound  is  a  significant  charge  reversal  on 
the  A1  site.  We  find  the  A1  atoms  to  carry  a  net  positive  charge  in  the  presence  of  hydrogen, 
whereas  they  carry  a  net  negative  charge  in  the  hydrogen-free  Ti3Al  alloy.  On  the  other 
hand,  independent  of  hydrogen  concentration,  the  net  charge  on  the  Ti  atom  is  always 
positive,  whereas  that  on  the  H  site  is  always  negative. 

The  calculated  total  charge  density  of  the  hexagonal  Ti12AlH6  cluster  [see  Fig.  7(a)] 
is  presented  in  Fig.  7(b).  These  calculations  confirm  the  inversion  of  the  net  charge  on 
the  A1  site,  depending  on  the  presence  or  absence  of  hydrogen  in  the  cluster.  The  total 
energy  difference  between  the  hydrogen  loaded  and  hydrogen-free  Tii2Al  cluster  indicate 
that  hydrogen  is  likely  to  be  absorbed  even  in  finite  Tii_xAlx  systems,  in  agreement  with 
experimental  observations  [8]. 

The  above  results  suggest  that  the  hydrogen-induced  decrease  of  the  shear  modulus 
in  Tix_xAlx  systems  may  be  correlated  with  the  net  charge  near  the  A1  site.  To  further 
elaborate  on  this  point,  we  compared  in  Figs.  9(a)  and  (b)  the  charge  distributions  in  the 
highest  occupied  molecular  orbitals  (HOMO,  corresponding  to  the  top  of  the  valence  band)  of 
Tii3H6  and  Tii2AlH6  clusters,  respectively.  The  geometry  and  the  relevant  shear  directions 
as  well  as  the  gliding  planes  in  these  systems  are  given  schematically  in  Fig.  9(c).  The 
calculated  HOMO  charge  distributions  suggest  the  presence  of  charge  lobes  perpendicular 
to  the  shear  direction  in  Tii3H6  clusters,  but  absence  of  such  bonding  charge  lobes  in  the 
Ti42AlH6  clusters.  The  absence  of  bonding  charge  regions  between  the  gliding  planes  may 
lie  at  the  origin  of  the  reduced  shear  modulus  in  Tii_xAlxHy  systems. 
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VI.  DISCUSSION 


As  already  mentioned  in  the  Introduction,  hydrogen-induced  embrittlement  occurring  in 
Tii_xAlx  alloys  is  a  very  complex  problem.  The  basic  questions  to  be  answered  are  (i)  what 
atomic-scale  processes  are  to  be  associated  with  embrittlement,  (ii)  which  changes  in  the 
electronic  structure  are  most  likely  to  trigger  these  processes,  and  (iii)  how  this  information 
could  be  used  to  reduce  or  inhibit  embrittlement  in  these  systems.  When  addressing  these 
problems,  we  have  to  bear  in  mind  that  the  response  of  the  system  to  applied  stress  depends 
crucially  both  on  the  static  structure  and  bonding  strength  in  the  system,  as  well  as  on  the 
deformation  dynamics  leading  to  mechanical  failure. 

The  first  basic  problem  to  be  addressed  is  the  equilibrium  structure  of  the  system  on  the 
atomic  scale.  Of  interest  is  not  only  the  equilibrium  geometry  of  an  ordered  monocrystalline 
alloy  of  a  given  composition,  but  also  the  possibility  of  phase  transitions  induced  by  internal 
and  external  stress,  phase  separation  in  the  bulk,  and  segregation  at  grain  boundaries  if  the 
system  is  polycrystalline.  Due  to  the  importance  of  defects  on  the  overall  stability  of  the 
system,  adequate  attention  must  be  paid  to  the  structure  of  dislocations  (including  partial 
and  superdislocations),  faulted  dipoles,  antiphase  boundaries  and  twin  boundaries,  and  the 
effect  of  other  defects  such  as  vacancies,  substitutional  and  interstitial  impurities  on  the 
equilibrium  structure. 

Even  more  important  than  the  structure  itself  is  the  nature,  strength,  and  anisotropy 
of  interatomic  bonding.  Of  primary  interest  in  the  Ti!_xAlxHy  systems  are  changes  in  the 
electronic  structure  and  the  changing  interplay  between  covalency  and  ionicity  of  interatomic 
bonding  as  a  function  of  composition,  in  particular  near  defect  sites.  Local  changes  of  the 
electronic  structure  will  affect  not  only  the  local  radial  compressibility  and  overall  strength 
of  interatomic  bonds,  which  is  reflected  in  the  bulk  modulus,  but  also  other  deformation 
modes  of  the  system,  such  as  those  reflected  in  the  shear  modulus.  Of  particular  importance 
for  the  overall  mechanical  stability  of  the  system  are  such  changes  near  defects  and  other 
weak  links. 
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Mechanical  deformations  under  applied  stress,  which  lead  to  failure,  proceed  via  inter¬ 
granular  and  transgranular  fracture  [36].  The  dynamics  of  the  transgranular  fracture,  which 
is  of  main  interest  in  this  work,  is  primarily  determined  by  the  ease  of  Peierls  relief  and  large- 
scale  dislocation  motion,  both  crucially  dependent  on  the  potential-energy  surfaces  and  the 
temperature-dependent  mobility  of  defects  (such  as  interstitial  hydrogen)  for  the  structures 
discussed  above.  Uncertainty  about  the  atomic-scale  processes  linked  with  embrittlement 
results  from  the  fact  that  whereas  brittle  fracture  is  typically  linked  to  the  inhibition  of 
dislocation  motion,  processes  such  as  the  Hydrogen  Enhanced  Local  Plasticity  of  a  stable 
system  can  effectively  masquerade  as  embrittlement  on  the  macroscopic  scale  [37-40].  This 
latter  process  will  be  discussed  in  more  detail  below. 

In  this  study,  we  attempted  to  answer  some  of  the  above  questions  by  calculating  from 
first  principles  the  equilibrium  geometry  and  stability,  the  electronic  structure,  and  the 
elastic  response  of  Tii_xAlxHy  systems.  We  used  three  different  implementations  of  the  ab 
initio  Density  Functional  formalism,  with  no  adjustable  parameters,  in  order  to  compare  the 
electronic  structure  and  bonding  strength  near  defects  to  those  of  similar  ordered  alloys  for 
a  wide  range  of  compositions.  The  computational  difficulties  for  these  systems  derived  from 
the  simultaneous  presence  of  Ti  atoms  with  strongly  localized  3 d  electrons,  A1  atoms  with 
much  more  delocalized  3 s  and  3p  electrons,  and  H  atoms  with  a  tendency  towards  covalent 
bonding,  as  well  as  the  large  charge  transfers  between  these  atoms.  We  found  all  our  results 
to  be  in  good  agreement  with  experimental  data  where  available. 

Our  calculations  as  well  as  an  inspection  of  the  phase  diagram  of  Tii_xAlx  alloys  indicate 
that  the  stability  deterioration  introduced  by  aluminum  and  hydrogen  atoms  are  not  caused 
by  local  structural  rearrangements  (“phase  transitions”),  but  must  originate  in  local  changes 
of  interatomic  bonding.  In  the  following,  we  will  discuss  which  changes  in  the  interatomic 
bonding  and  the  underlying  electronic  structure  were  found  to  correlate  with  the  stability 
reduction  occurring  in  Tii_xAlx  alloys  upon  hydrogen  loading.  Even  though  such  static 
considerations  have  been  proven  inadequate  to  describe  quantitatively  the  stability  limits  of 
perfect  crystals  at  finite  strain  [41],  we  use  them  as  a  first  step  towards  the  understanding 
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of  the  alleged  “hydrogen  embrittlement”  in  these  systems. 

The  first  question  to  be  answered  is,  whether  the  stability  of  metallic  Ti  is  significantly 
reduced  by  the  simultaneous  presence  of  A1  and  H  atoms.  Our  results  for  the  cohesive  en¬ 
ergies  Ecoh  (with  respect  to  isolated  atoms),  presented  in  Table  I,  indicate  that  presence  of 
hydrogen  stabilizes  all  the  metallic  systems.  This  is  true  when  comparing  ^^(TiH)  and 
Ecoh( TiH2)  to  .Eco/^Ti),  Ecoh(AlK)  to  J5coh( Al),  and  Eco/l(Ti3AlH4)  to  (Ti3Al).  Among 
these  systems,  we  found  the  allegedly  brittle  Ti3AlH4  compound  to  be  most  stable,  and 
hence  brittle  behavior  to  anticorrelate  with  increasing  stability.  Our  results  for  the  for¬ 
mation  energies  Ef^m  (with  respect  to  the  elements  in  the  standard  state),  presented  in 
Table  II,  indicate  that  also  the  presence  of  Al  in  Ti  leads  to  a  stable  Tii_xAlx  alloy,  even 
though  hydrogen  would  prefer  to  be  close  to  the  Ti  sites  in  systems  with  low  hydrogen  con¬ 
centrations.  The  Mulliken  population  analysis  in  Table  III  indicates  that  the  large  stability 
of  the  Tii_xAlxHy  systems  is  a  consequence  of  the  large  charge  transfer  between  the  sites, 
resulting  in  a  partly  ionic  bonding. 

After  discarding  a  possible  link  between  embrittlement  and  a  reduced  stability,  we  in¬ 
vestigated  the  correlation  between  reduced  mechanical  strength  and  the  elastic  response  of 
the  different  Tix_xAlxHy  systems.  Our  results  for  the  bulk  modulus  are  also  summarized 
in  Table  I.  These  results  confirm  that  with  its  bulk  modulus  being  more  than  twice  that 
of  its  pure  constituents,  Ti3Al  is  indeed  much  stiffer  and  stronger  than  pure  bulk  Ti  or  Al. 
Our  results  indicate  that  the  presence  of  hydrogen  has  only  a  negligible  effect  on  the  bulk 
modulus  of  all  Tii_xAlx  systems  investigated.  Hence  we  conclude  that  the  embrittlement  in 
Tii_xAlxHy  systems  is  not  reflected  in  changes  of  the  bulk  modulus. 

Finally,  we  investigated  the  elastic  response  of  Tii_xAlxHy  systems  with  respect  to  ap¬ 
plied  shear  stress,  noting  that  a  decreasing  shear  modulus  C44  would  likely  enhance  the 
creation  and  mobility  of  dislocations  and  hence  lead  to  a  locally  enhanced  plasticity.  We 
calculated  C44  using  total  energy  differences,  and  verified  the  linearity  of  the  response  for 
small  deformations.  Due  to  the  large  size  of  the  unit  cells  (which  are  implicitly  assumed  to  be 
orthogonal  in  our  bulk  LDA  code)  in  distorted  single  crystals,  we  performed  the  calculations 
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of  C44  on  small  Tii_xAlxHy  clusters  with  an  hep  geometry.  There,  we  were  not  restricted  to 
periodic  structures,  and  could  consider  even  very  small  deformations  along  the  a;- direction 
perpendicular  to  the  c  axis.  A  good  agreement  between  cluster-based  and  bulk-based  cal¬ 
culations,  verified  in  the  case  of  Al(fcc),  indicates  that  our  results  axe  not  dominated  by 
surface  effects. 

Our  results  for  the  elastic  constant  C44  in  Tii_xAlxHy  systems,  for  a  shear  of  the  (0001) 
planes  in  response  to  a  shear  stress  along  the  [1000]  direction,  are  summarized  in  Table  I. 
As  in  the  case  of  the  bulk  modulus,  we  find  that  also  the  shear  modulus  of  Ti  increases 
significantly  upon  alloying  with  Al.  In  strong  contrast  to  the  bulk  modulus,  however,  we 
find  a  strong  decrease  of  the  shear  modulus  in  Tii_xAlx  alloys  upon  hydrogen  loading,  namely 
from  C,44(Ti3Al)=  101  GPa  to  C44(Ti3AlH4)=  67  GPa.  We  conclude  that  the  decrease  of  the 
shear  modulus  C44  of  Ti,  which  occurs  only  while  hydrogen  and  aluminum  are  simultaneously 
present  in  the  metal  matrix,  is  closely  related  to  the  reduced  mechanical  strength  of  Tii_xAlx 
systems  upon  hydrogen  loading.  This  interpretation  is  also  consistent  with  our  result  of 
Table  I  that  in  absence  of  aluminum,  the  introduction  of  H  into  the  Ti  matrix  leads  to 
an  increase  of  the  shear  modulus  from  C44(Ti)=  43  GPa  to  C,44(TiH)=  52  GPa,  in  close 
agreement  with  the  analogous  behavior  observed  in  NbHy  systems  [42]. 

Less  easy  to  interpret  are  our  results  for  the  elastic  response  to  applied  shear  stress 
along  the  [1120]  direction,  which  results  in  a  shear  of  the  (1100)  planes.  Even  though 
these  planes  are  known  to  form  the  preferential  slip  system  in  the  Tii_xAlx  alloys  [43],  our 
preliminary  results  indicate  relatively  large  values  of  C4 4,  ranging  from  (744(Ti)=  137  GPa, 
to  (744(Ti3Al)=  129  GPa,  and  C,44(Ti3AlH4)=  143  GPa. 

The  hydrogen-induced  reduction  of  the  shear  modulus  along  the  [1000]  direction  in 
Tii_xAlx  alloys,  and  the  anticipated  increased  mobility  of  dislocations  in  these  systems,  are 
consistent  with  the  “Hydrogen  Enhanced  Local  Plasticity”  (HELP)  mechanism  of  failure. 
This  microscopic  mechanism  had  originally  been  proposed  by  C.  Beachem  [37]  to  explain 
fracture  of  hydrogen-loaded  steel  under  tensile  stress,  and  has  been  more  recently  used  to 
explain  a  similar  behavior  in  other  high-strength  Al-based  alloys  [39].  In  systems  under 
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tensile  stress,  defects  axe  often  present  near  strain  singular  points  such  as  crack  tips.  The 
initial  growth  and  propagation  of  these  cracks  is  accelerated  by  a  rapid  local  accumulation  of 
hydrogen,  which  leads  to  a  local  ductility  enhancement  near  these  crack  tips.  Under  tensile 
stress,  fracture  propagation  is  accelerated  at  such  hydrogen  saturated  weak  links,  while  the 
bulk  of  the  matrix  is  essentially  unperturbed.  This  mechanism  leads  to  a  fracture  morphol¬ 
ogy  which  is  indistinguishable  from  brittle  fracture,  since  brittle  fracture  implies  no  crystal 
shape  changes  outside  the  crack  region,  whereas  a  typical  ductile  fracture  is  accompanied 
by  large-scale  shape  deformations. 

The  next  important  problem  is  to  determine  the  electronic  origin  of  the  strong  reduction 
of  the  shear  modulus  in  Tii_xAlx  alloys  upon  hydrogen  loading,  which  we  linked  to  the 
“hydrogen  embrittlement”  of  these  systems.  In  particular,  we  want  to  identify  the  electronic 
properties  which  would  reflect  changes  of  C44. 

One  of  the  most  suitable  quantities  to  reflect  bonding  changes  in  a  given  system  is  the 
electronic  density  of  states.  Indeed,  we  find  the  electronic  density  of  states  of  Ti3Al,  shown 
in  Fig.  8(a),  to  differ  significantly  from  that  of  Ti3AlH4,  given  in  Fig.  8(b).  Even  though 
the  increased  population  of  s-states,  apparent  from  the  comparison  of  Figs.  8(a)  and  (b), 
may  be  indicative  of  a  decreased  bonding  anisotropy  which  would  cause  a  lowering  of  C44, 
a  definitive  conclusion  along  these  lines  is  difficult  to  make  based  on  these  results. 

The  difficulty  to  extract  information  about  the  shear  modulus  from  the  density  of  states 
results  form  the  fact  that  bonding  anisotropy  is  only  very  indirectly  reflected  in  this  quan¬ 
tity.  From  this  point  of  view,  an  observable  which  should  reflect  more  directly  the  bonding 
anisotropy  is  the  spatial  distribution  of  electrons.  We  would  expect  the  most  pronounced 
changes  due  to  the  presence  of  hydrogen  to  occur  near  a  substitutional  A1  impurity  in 
bulk  Ti(hcp).  Our  corresponding  results  for  a  T^Al  and  a  Tii2AlH6  clusters  are  shown  in 
Figs.  3(b)  and  7(b),  respectively.  These  results  indicate  the  predominant  localization  of  the 
electronic  charge  on  the  atomic  sites.  Even  the  more  detailed  charge  density  contour  plot 
in  Fig  3(c)  for  the  hydrogen-free  system  does  not  show  any  significant  charge  accumulation 
in  the  bond  regions,  or  a  bonding  anisotropy  of  the  charge  density.  Hence  we  find  no  simple 
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indication  of  a  reduced  shear  modulus  in  the  total  charge  density  distribution. 

A  first  indication  of  significant  electronic  structure  differences  between  the  hydrogen-free 
Ti3Al  alloy,  which  has  a  large  shear  modulus,  and  the  hydrogen-loaded  Ti3Al  alloy,  which 
has  a  low  shear  modulus,  was  observed  in  the  Mulliken  populations  of  the  sites.  Our  results, 
presented  in  Table  III,  show  that  in  all  the  systems  investigated  in  this  study,  hydrogen 
always  carries  a  small  negative  excess  charge,  while  Ti  always  carries  a  small  positive  net 
charge.  Of  particular  interest,  however,  is  a  sign  reversal  of  the  net  charge  on  the  A1 
site,  depending  on  whether  or  not  hydrogen  is  present  in  the  metal  matrix.  We  can  easily 
imagine  that  a  partial  ionic  bonding  between  oppositely  charged  A1  and  Ti  sites  may  provide 
additional  stabilization  within  the  hydrogen-free  system  with  respect  to  shear  distortions. 
This  effect  is  absent  in  the  hydrogen-loaded  Ti3Al  alloy  where  all  metal  sites  carry  a  net 
positive  charge.  Hence  on  the  electronic  level,  we  may  tentatively  link  the  same  sign  of 
the  excess  charge  on  the  metal  impurity  and  its  neighbors  in  the  matrix  with  a  propensity 
towards  lowering  the  shear  modulus. 

Perhaps  more  instructive  that  the  net  charge  on  the  individual  sites,  or  the  distribution 
of  the  total  charge  density,  should  be  an  inspection  of  the  valence  charge  distribution.  In 
clusters,  as  well  as  in  molecules,  it  is  the  charge  distribution  in  the  Highest  Occupied  Molec¬ 
ular  Orbital  (HOMO)  which  reflects  the  bonding  nature,  the  stability,  and  particularly  weak 
links  in  the  system.  The  HOMO  charge  distribution  in  the  hexagonal  clusters  Tii3H6  and 
Tii2AlH6  is  shown  in  Figs.  9(a)  and  (b),  respectively.  The  HOMO  of  the  Al-free  system 
displays  lobes  of  charge  in  the  direction  of  the  c-axis  which  connect  as  “bridges”  the  (0001) 
planes  in  the  structure.  Such  valence  charge  “bridges”  connecting  the  (0001)  planes  are 
even  stronger  in  the  hydrogen-free  ordered  alloy  Ti3Al  with  a  large  value  of  C44,  and  near 
substitutional  A1  impurities  in  Ti,  represented  by  a  Ti^Al  cluster.  On  the  other  hand,  as 
shown  in  Fig.  9(b),  these  lobes  or  “bridges”  are  absent  near  substitutional  A1  impurities 
in  the  hydrogen-loaded  Ti  matrix.  Hence  in  the  particular  case  of  Tii_xAlxHy  alloys,  we 
interpret  a  connected  three-dimensional  mesh  of  charge  density  in  the  HOMO,  such  as  that 
found  for  Tii3H6  and  Tii2Al,  as  an  indication  for  a  shear-resistant  structure.  On  the  other 
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hand,  presence  of  charge-depleted  (0001)  planes  in  the  HOMO  charge  distribution,  such  as 
in  the  case  of  Ti^AlHe,  can  be  plausibly  connected  with  a  lowered  resistance  towards  shear 
along  a  direction  perpendicular  to  the  c-axis. 

VII.  CONCLUSIONS 

We  used  a  combination  of  different  ab  initio  computational  methods  in  an  attempt  to 
understand  the  adverse  effect  of  hydrogen  on  the  mechanical  stability  of  Tii_xAlx  alloys, 
also  known  as  “hydrogen  embrittlement”.  Our  studies  focussed  on  the  equilibrium  structure 
and  the  elastic  constants  of  bulk  Ti,  Al,  the  ordered  alloy  Ti3Al,  as  well  as  their  hydrides. 
Calculations  for  ordered  systems  were  performed  using  a  bulk  method  with  local  orbitals  and 
ab  initio  pseudopotentials.  Results  for  the  electronic  structure  of  Ti  near  atomic  vacancies, 
substitutional  Al  and  interstitial  H  impurities  were  obtained  using  the  cluster  approach  and 
the  Linear  Muffin- Tin  Orbitals  in  the  Atomic  Spheres  Approximation  (LMTO-ASA)-Green 
function  formalism. 

The  complementary  results  of  these  calculations  indicate  that  the  reduced  stability  of 
Tii_xAlx  alloys  is  reflected  in  a  significant  reduction  of  the  shear  modulus  C44  of  these 
systems  upon  hydrogen  loading.  On  the  electronic  level,  we  found  the  decrease  of  the  shear 
modulus  to  be  reflected  in  the  distribution  of  the  electrons  at  the  top  of  the  valence  band, 
or  in  the  HOMO  orbital  in  the  corresponding  clusters.  While  this  valence  charge  forms  a 
relatively  rigid  mesh  in  Ti3Al  and  TiH,  we  observe  a  valence  charge  depletion  in  the  (0001) 
planes  of  the  lattice  in  the  case  that  both  Al  and  H  are  present  in  the  Ti  matrix.  It  appears 
as  plausible  that  this  valence  charge  depletion  may  reduce  the  resistance  of  hydrogen-loaded 
Tii_xAlx  alloys  towards  shear  perpendicular  to  the  c-axis.  Consequently,  these  planes  could 
act  as  preferential  gliding  planes  in  perfect  crystals  or  preferential  planes  for  dislocation 
motion  in  imperfect  crystals,  resulting  in  a  plastic  deformation  of  the  titanium  matrix  near 
aluminum  impurities  surrounded  by  hydrogen.  This  interpretation  would  be  consistent 
with  the  Hydrogen-Enhanced  Local  Plasticity  Mechanism  discussed  above.  We  note  that 
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all  our  theoretical  findings,  as  well  as  their  above  interpretation,  are  consistent  with  the 
experimental  observation  that  hydrogen  embrittlement  occurs  only  in  presence  of  A1  in  the 
Ti  matrix  [10]. 

Of  course,  these  static  results,  obtained  within  the  linear  response  regime,  are  far  from 
describing  the  complex  nature  of  fracture  in  realistic  systems.  Nevertheless,  they  provide 
us  with  a  first  understanding  of  the  atomic-level  processes  which  are  likely  to  occur  in 
hydrogen-loaded  Tii_xAlx  systems  subjected  to  tensile  stress,  and  will  form  a  basis  for  future 
dynamical  calculations.  Such  studies  are  likely  to  address  the  effect  of  extended  stacking 
faults  and  microvoids  on  the  overall  stability,  the  dynamics  of  dislocation  motion,  as  well 
as  the  mechanism  of  crack  propagation  in  these  alloys.  We  believe  that  a  fundamental 
understanding  of  atomic-  and  electronic-level  processes  linked  with  the  adverse  effect  of 
hydrogen  on  the  mechanical  stability  of  Ti!_xAlx  systems  could  prove  to  be  very  useful 
in  the  exploration  of  materials  treatment  techniques  and  additives  which  would  inhibit  or 
suppress  hydrogen  embrittlement  in  these  important  alloys. 
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TABLES 


TABLE  I.  Calculated  and  observed  ground  state  properties  of  bulk  Al,  Ti,  Ti3Al,  and  their 
hydrides:  lattice  constant  a,  cohesive  energy  (with  respect  to  isolated  atoms)  Ecoh,  bulk  modulus 


B ,  and  shear  modulus  <744. 


Al 

Ti(fcc) 

Ti(hcp) 

Ti3Al 

A1H 

TiH 

TiH2 

Ti3AlH4 

a  (a.u.) 

LDAa 

7.80 

7.73 

5.55 

7.76 

8.25 

8.19 

8.63 

8.14 

Expt. 

7.66 

7.73 

5.58 

7.68 

8.41 

8.22 

E^  (eV) 

LDA 

-3.51 

-3.70 

-3.70 

-15.15 

-4.99 

-6.15 

-8.81 

-24.19 

Expt. 

-3.39 

-4.70 

-4.70 

B  (GPa) 

LDA 

71 

94 

100 

280 

70 

72 

74 

290 

Expt. 

73 

106 

106 

C44  (GPa) 

LDAb 

37 

43 

101 

15 

52 

67 

Expt. 

32 

46 

a  Results  obtained  using  the  LMTO-ASA  method  are  discussed  in  footnote  [24]. 
b  Calculations  are  for  shear  within  the  (0001)  plane,  along  the  [1000]  direction  [see  Fig.  9(c)]. 


TABLE  II.  Formation  energies  Efo,.m  (eV)  [29,31]  of  Tii xAlxHy  systems. 


Ti3Al 

A1H 

TiH 

TiH2 

Ti3Affl4 

-0.54 

0.87 

-0.10 

-0.40 

0.45 

TABLE  m.  Mulliken  excess  charges  (e  )  on  individual  sites  in  Tii xAlxHy  systems. 


A1H 

TiH 

TiH2 

Ti3Al 

Ti3AlH4 

Ti 

0.28 

0.08 

0.30 

0.65 

Al 

0.70 

-0.89 

0.43 

H 

-0.70 

-0.28 

-0.04 

-0.59 
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FIGURES 


FIG.  1.  Electronic  density  of  states  of  bulk  (a)  Al(fcc),  (b)  Ti(hcp),  and  (c)  Ti(fcc),  calculated 
using  ab  initio  pseudopotentials  and  a  local  orbital  basis.  The  bold  solid  line  gives  the  total  density 
of  states,  the  solid  line  the  partial  density  of  states  (PDOS)  of  s-electrons,  the  dashed  line  the  PDOS 
of  p-electrons,  and  the  dotted  line  the  PDOS  of  d-electrons.  Note  the  different  vertical  scales  used 
in  (a),  (b)  and  (c). 

FIG.  2.  Electronic  density  of  states  of  the  ordered  bulk  alloy  Ti3Al  in  the  cubic  phase,  calculated 
using  ab  initio  pseudopotentials  and  a  local  orbital  basis,  (a)  Total  density  of  states  per  Ti3Al 
formula  unit,  (b)  the  local  density  of  states  (LDOS)  at  the  A1  site,  and  (c)  the  LDOS  at  the  Ti 
site.  Note  the  different  vertical  scales  used  in  (a),  (b)  and  (c). 

FIG.  3.  Geometry  and  charge  distribution  in  a  very  dilute  TiAl  alloy,  represented  by  Ti(hcp) 
with  a  substitutional  A1  impurity,  (a)  Geometry  of  the  Ti12Al  core  cluster,  (b)  3-D  representation  of 
the  total  charge  density  surface  in  the  Ti^Al  cluster,  corresponding  to  p  =  0.1  e-/a.u.3,  calculated 
using  the  DMol  all-electron  cluster  method,  (c)  Results  of  the  same  calculation,  presented  as  a  2-D 
contour  plot  in  the  (0001)  plane,  which  is  shown  by  the  the  dashed  line  in  (a)  and  contains  the  A1 
impurity,  (d)  Total  charge  density  in  the  dilute  TiAl  alloy,  calculated  using  the  LMTO-ASA-Green 
function  method,  and  represented  in  the  same  manner  and  geometry  as  in  (c).  The  ratio  of  two 
consecutive  contours  in  (c)  and  (d)  is  p(n  +  l)/p(n)  =  1.5. 

FIG.  4.  Electronic  structure  of  a  very  dilute  TiAl  alloy,  represented  by  Ti(hcp)  with  a  substi¬ 
tutional  A1  impurity,  calculated  using  the  LMTO-ASA-Green  function  method.  The  local  density 
of  states  at  (a)  the  A1  site  and  (b)  the  neighboring  Ti  site  is  given  by  the  solid  line. 

FIG.  5.  Geometry  and  charge  distribution  in  TiH,  with  hydrogen  atoms  occupying  the  octahe¬ 
dral  interstitial  sites  of  the  hep  lattice,  (a)  Geometry  of  a  Tii3H6  core  cluster,  (b)  3-D  represen¬ 
tation  of  the  total  charge  density  surface  in  the  Tii3H6  cluster,  corresponding  to  p  =  0.1  e~/a.u.3, 
calculated  using  the  DMol  all-electron  cluster  method. 
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FIG.  6.  Comparison  of  electronic  densities  of  states  in  bulk  Ti  and  TiH2.  (a)  Density  of  states 
of  Ti(fcc),  (b)  density  of  states  of  TiH2(fcc),  and  (c)  local  density  of  states  on  the  Ti  site  in 
TiH2(fcc),  calculated  using  ab  initio  pseudopotentials  and  a  local  orbital  basis.  The  bold  solid  line 
gives  the  total  density  of  states,  the  solid  line  the  partial  density  of  states  (PDOS)  of  s-electrons, 
the  dashed  line  the  PDOS  of  p-electrons,  and  the  dotted  line  the  PDOS  of  d-electrons. 

FIG.  7.  Geometry  and  charge  distribution  in  the  hydride  of  a  very  dilute  TiAl  alloy,  rep¬ 
resented  by  a  Ti(hcp)  matrix  with  a  substitutional  A1  impurity  and  hydrogen  atoms  occupying 
the  octahedral  interstitial  sites  of  the  hep  lattice,  (a)  Geometry  of  a  Ti^AlHe  core  cluster,  (b) 
3-D  representation  of  the  total  charge  density  surface  in  the  Ti^AlHe  cluster,  corresponding  to 
p  =  0.1  e~/a.u.3,  calculated  using  the  DMol  all-electron  cluster  method. 

FIG.  8.  Comparison  of  the  total  electronic  densities  of  states  in  bulk  (a)  TisAl(fcc)  and  (b) 
Ti3AlH4,  calculated  using  ab  initio  pseudopotentials  and  a  local  orbital  basis.  The  bold  solid  line 
gives  the  total  density  of  states,  the  solid  line  the  partial  density  of  states  (PDOS)  of  s-electrons, 
the  dashed  line  the  PDOS  of  p-electrons,  and  the  dotted  line  the  PDOS  of  d-electrons. 

FIG.  9.  Charge  distribution  in  the  highest  occupied  molecular  orbitals  (HOMO)  of  (a)  Ti^He 
and  (b)  Ti^AlHe  clusters,  both  based  on  the  bulk  Ti(hcp)  structure,  calculated  using  the  DMol 
all-electron  cluster  method.  A  3-D  representation  of  the  total  charge  density  surface  in  (a)  and 
(b)  is  shown  at  the  value  p =  0.001  e-/a.u.3.  The  geometry  has  been  given  in  Figs.  5(a)  and  7(a), 
respectively.  The  arrows  indicate  possible  applied  shear  stress,  (c)  Schematic  drawing  of  the  shear 
directions  and  gliding  planes  in  Tii_xAlx  systems. 
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